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ON  THE  CONVECTIVE  PROPERTIES  OF  THE  RESISTIVE  HOSE  INSTABILITY 
IN  AN  INTENSE  RELATIVISTIC  ELECTRON  BEAM 


I.  Introduction 

The  resistive  hose  instability  is  a  dominant  instability  of  an  intense 

relativistic  electron  beam  propagating  in  a  weakly  ionized  or  initially 

1  2 

unionized  gas.  It  has  been  subject  to  a  variety  of  experimental  *  and 
3—8 

theoretical  studies  in  recent  years.  Although  the  formulation  of  the 
physics  involved  in  beam  propagation  is  a  difficult  problem,  even  as  regards 
the  construction  of  an  equilibrium,  remarkable  progress  has  been  made  in 
understanding  the  hose  instability  by  means  of  highly  simplified  models.  At 
present,  there  is  substantial  qualitative  correlation  among  simple  analytic 
theory,  experiment,  particle  simulation,  and  other  techniques  of  numerical 
analysis . 

There  are  certain  minimum  requirements  that  must  be  satisfied  by  any 

acceptable  model  of  the  hose  instability.  For  example,  the  early  work  based 

9  10 

on  the  "rigid  beam  model"  or  on  the  assumption  of  a  constant  betatron 

frequency  for  all  beam  electrons  indicated  that  the  instability  is  absolute  in 

the  frame  co-moving  with  the  beam.  More  realistic  models,  such  as  the  spread 

mass  model, J  or  the  multicomponent  model,  indicate  that  the  Instability  is  of 

a  convective  nature.  However,  direct  examination  of  the  dispersion  relation 

does  not  yield  much  information  on  the  minimum  velocity  of  convection  of  the 

amplifying  waves,  which  may  be  important  if  one  wants  to  assess  the  effect  of 

the  beam  head  erosion  on  the  evolution  of  the  instability.  Furthermore, 

recent  work4* 3,0  has  emphasized  that  a  plasma  return  current  destabilizes  the 

hose  instability,  and  even  leads  to  nonvanishing  spatial  amplification  rates 

at  very  low  frequency,  thereby  leading  to  the  speculation  that  the  hose 

instability  may  indeed  become  absolute  for  sufficiently  large  return  current 

fraction. 
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In  the  present  paper,  we  study  the  space-time  evolution  of  the  resistive 

hose  instability.  Ue  include  the  effects  of  plasma  return  current  as  well  as 
a  self-consistent  treatment  of  the  plasma  electrical  conductivity  increase 
along  the  beam  axis,  due  (for  example)  to  additional  beam-collisional 
ionization  of  the  gas.  The  analysis  is  essentially  one-dimensional  (axial)  in 

Q 

space,  and  is  thus  restricted  mainly  to  the  beam  body,  the  region  behind  the 
pinch  point  >  where  the  beam  current  and  radius  may  be  assumed  to  be 
roughly  coastant,  and  dependence  on  radial  position  is  easily  separated  out. 
Under  the  simplifying  assumptions  enumerated  at  the  beginning  of  Section  II, 
we  show  that  the  resistive  hose  instability  is  always  convective,  regardless 
of  the  level  of  the  return  current.  This  conclusion  results  from  our 
examination  of  the  Green's  function  J  constructed  from  the  dispersion  relation 

O 

which  was  derived  by  Lampe,  Sharp,  and  Hubbard0  for  the  spread  mass  model.  In 
Section  III,  we  extend  our  analysis  to  calculate  in  more  detail  the  evolution 
of  amplifying  disturbances  based  on  an  ad  hoc  simplification  of  the  dispersion 
relation.  We  find  that  there  is  a  minimum  speed  of  propagation  for  amplifying 
waves.  This  minimum  speed  is  model  dependent,  but  is  typically  at  most  half 
of  the  speed  associated  with  the  maximum  convective  growth.  This  last 
statement  holds  throughout  the  vast  parameter  space  surveyed.  The  maximum 
spatial  amplification  rate  is  just  the  value  expected  from  the  dispersion 


relationship 


II.  The  Model 


In  this  section,  we  describe  our  model  of  the  relativistic  electron 
beam.  We  show  here  that  the  resistive  hose  instability  is  always  convective, 
when  the  spread  mass  model  is  adopted.  We  postpone  to  the  next  section  a 
calculation  of  the  propagation  characteristics. 

Our  analysis  applies  to  the  main  body  of  the  beam,  at  a  distance 
sufficiently  far  from  the  beam  head.  The  following  simplifying  assumptions 
are  introduced:  (a)  The  electron  beam  is  highly  relativistic  and  paraxial, 
with  long  betatron  wavelength  compared  with  the  beam  radius.  (b)  The 
electrical  conductivity  of  the  plasma  channel  is  high  enough  that  the  charge 
relaxation  length  is  short  compared  with  beam  radius.  (c)  The  beam  current 
density  Jj,0(r)»  P^asma  current  density  Jp0(r)»  and  the  conductivity  0Q(r),  all 
assume  the  same  Bennett  profile  and  are  constant  or  slowly  varying  in  z  and 
t.8  (d)  The  plasma  return  current  Ip  is  a  constant  fraction  f  of  the  beeam 

current  Ib,  with  the  convention  that  f  is  negative  when  the  plasma  current 
flows  in  the  opposite  direction  to  the  beam  current. ^»^»8  (e)  The  rate  of 

increase  of  the  equilibrium  electrical  conductivity  0  q  is  proportional  to  the 
local  current  density,8’^8  due  to  beam-impact  ionization  of  the  gas. 
Mathematically,  the  last  assumption  can  be  written  AaQ/Ax,  m  k  J^,  where  5  is 
the  axial  distance  measured  from  the  beam  head  toward  the  beam  tail  and  k  is 

£  O 

assumed  to  be  constant.  An  important  consequence  *  of  this  assumption  is 
that  (in  ?),  rather  than  5,  is  an  ignorable  coordinate  for  the  governing 
equations  of  the  equilibrium  quantities.  In  addition,  (f)  perturbations  of 
the  electrical  conductivity  due  to  the  effects  of  hose-perturbed  beam  currents 
are  included  self-consistently ,  as  described  in  Ref.  8.  Assumptions  (a)  -  (f) 
have  been  adopted  in  the  literature  cited. 

Instead  of  using  the  independent  space-time  variables  (z,t)  in  the 


laboratory  frame,  with  z  being  the  direction  of  beam  propagation,  we  shall  use 
(z,?)  as  Independent  variables  In  the  description  of  the  hose  instability. 
Here,  (z,;)  is  related  to  (z,t)  through  the  transformation 


5  -  vfat  -  z,  (1) 

Z  -  z,  (2) 

and  ;  may  be  regarded  as  the  distance  measured  from  the  beam  head  toward  the 

beam  body,  and  is  thus  a  "space-like"  variable,  while  z  becomes  a  "time-like" 

variable  for  a  fixed  distance  5  away  from  the  beam  head.  In  Eq.  (1),  v^  is 

the  beam  velocity  and  is  very  close  to  c,  the  speed  of  light.  As  mentioned  in 

the  last  paragraph,  in  5  is  an  ignorable  coordinate.  Small  amplitude 

perturbed  quantities  are  then  expected  to  vary  as  exp  [-  i  ^roz  -  1  u>  in  5] 

where  nrQ  is  the  "frequency"  associated  with  "temporal"  variation  in  z  and  to 

is  the  "wavenumber"  associated  with  "spatial"  variation  in  in  5.  The 

dispersion  relation  for  the  hose  instability  then  relates  the  frequency 

to  the  wavenumber  u*  Information  concerning  the  absolute/convective  behavior 

of  the  instability,  as  well  as  details  of  pulse  propagation  characteristics, 

13 

may  be  extracted  from  such  a  dispersion  relationship. 

It  is  convenient  to  use  the  normalized  independent  variables  Z,  ip , 
defined  by 


Z 


So2 


(3) 


<I>  -  (1/A)  in  (?/ a) 


(4) 


VJV.V. 


where  a  is  the  radius  of  the  Bennett  profile,  o  =  ( 2 1 1  !/I  is  the 

betatron  wavenumber,  and  X  =  <  Ib/2c  is  a  dimensionless,  positive  constant. 

For  propagation  in  air,  a  reasonable  estimate®  is  X  -  0.044  1^  with  Ib  in 
kilo-amperes.  Here,  the  net  current  In  is  the  sum  of  the  beam  current  and 
the  plasma  current  Ip,  and  1^  5  (17  g  y)  kA  is  the  Alfven-Lawson  current.  The 
normalized  "frequency"  n  =  ^rQ/0go  and  the  normalized  "wavenumber"  Xm  may  then 
be  introduced  to  the  linearized  quantities,  all  of  which  behave  according  to 

Linearized  quantity  ~  e  *  ^  ^  5) 

Under  the  assumptions  stated  at  the  beginning  of  this  section,  the  spread 
mass  model  yields  the  dispersion  relationship  D(  *  0,  where 

D(SJ.uj)  =»  [o  -  uij  (0)  ]  [to  -  u>2 (Q)  ]  (6) 

and  u»^,  (jj2  are  defined  as 

X“l,2  "  '  1  F  \  [R(n)]1/2.  (7) 

In  Eq.  (7), 


F(fl) 


=  X  + 


f  -  G(n) 
1  +  f 


(8) 


and 


H(n)  =  f2(0)  +  4  xG(a), 


(9) 


OTKi 


G(n)  =  /  dn 


~  T|  /  it 


n  -  n 


-  6Q2  {(i  -  n2)  +  fl2(  1  -  fl2)  in  (10) 

Q 

3 

Is  the  spread  mass  function  introduced  by  Lee.  The  derivation  of  the 

dispersion  relationship  (6)  is  given  by  Lampe,  et  al.  in  Ref.  8.  Before  we 

calculate  the  spatial  amplification  rate  Im(Xa))  as  a  function  of  the 
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frequency  we  shall  show  that  the  hose  instability  described  by  the 
dispersion  relationship  (6)  is  convective,  regardless  of  the  fraction  of  the 
return  current  f. 

The  following  properties  regarding  the  dispersion  relationship  will  be 
needed:  The  functions  G(n)  and  F(fl),  are  analytic  for  all  Q  with  Imfi  >  0. 

The  branch  points  of  G(n)  and  F(Q)  are  located  at  f|  «  0,  fi  *  +  1.  Branch  cuts 
originating  from  these  branch  points  are  introduced  to  render  these  functions 
analytic  in  the  upper  half  n  plane.  Futhermore,  |G(n)l  and  |F(fi)|  are  bounded 


as  Ifil  +  «. 


To  show  that  the  resistive  hose  instability  as  described  by  Eq.  (6)  is  a 
convective  instability,  we  need  only  to  consider  the  response  function 


R(Z,<|>) 


/  dn  s(n)  e_iaz 

L 


f 


-  iXaup 

d"  i - - 

D(n,u) 


(U) 


where  S(ft)  denotes  the  spectrum  of  a  source  of  excitation  and  L  is  the 
Bromwich  contour  on  which  Im(fi)  is  sufficiently  large  and  positive.  Absence 
of  absolute  instabi  'ty  is  '  .ablished  if  it  can  be  shown  that  R(Z,  <j>)  -»■  0  as 
"time"  Z  +  »,  at  any  '  •  xed  value  of  "position”  \j>,  for  any  well-behaved 


function  S(q).  J 

The  u)  -  integral  in  (11),  designated  as  M(q),  can  be  carried  out 
explicitly  by  the  method  of  residues: 


M(q) 


! 


dZ  e~ 

[oi  -  ^(q) ]  [ti>  -  1 


(12) 


It  can  easily  be  verified  that  as  Imq  ■*■  »,  >  -  i  and  +  -  iX.  Thus,  (12) 

yields 


M(a)  - 


{ 


0 

4iriexp[-  |  F(n)] 


sinh[|  H1/2(q) ] 


H1/2(«) 


,  for  <  0 
,  for  i|>  >  0. 


(13) 


Equation  (13)  is  valid  for  all  q  with  Imq  +•  «.  However,  M(q)  as  given  by  (13) 

may  be  analytically  continued  for  all  q  with  Imq  >  0  by  observing  that  M(q) 

is  singular  only  at  the  branch  points  of  F(q),  viz,  q  -  0,  +1,  and  at  the 

roots  q  of  H.  The  singularities  at  ft  are  removable  singularities, 
s  s 

Thus,  M(q)  is  analytic  for  all  Imq  >  0,  except  at  the  branch  points 
q  **  0,  +  1.  The  inverse  Laplace  transform  integral  (11)  may  then  be  evaluated 
along  the  real  q  axis: 


R(Z ,\j>)  -  /  dqS(q)e"inZ  M(q). 


(14) 


Now,  for  any  fixed  value  of  ij»,  !M(q)|  is  bounded  for  all  real  values  of  q. 
Thus,  if  the  source  function  S(q)  satisfies  certain  mild  restrictions,  (14) 


yields  R(Z » \p)  +  0  as  !  +  »  by  the  Rieraann-Lebesgue  Lemma.  Note  that  this 
result  is  independent  of  f,  the  fraction  of  the  plasma  return  current.  We 
have  thus  proven  that  the  hose  instability  is  always  convective. 

We  close  this  section  with  the  following  remarks,  (a)  The  disturbance 
always  propagates  from  the  beam  head  toward  the  beam  tail.  This  obvious  fact 
is  clearly  indicated  by  Eq.(13),  which  shows  that  no  perturbation  on  the  beam 
can  move  faster  than  the  beam  head  in  the  laboratory  frame.  (b)  Convective 
growth  is  present,  for  if  we  let  1)1  *  VZ  in  Eq.  (13b)  for  some  fixed  velocity 
V,  M(fi)  may  become  exponentially  large  as  Z  +  ».  We  shall  examine  this  point 
in  more  detail  in  the  next  section.  (c)  The  detailed  study  of  the  evolution 
of  the  response  function  R  as  given  by  Eq.  (14)  is  quite  involved.  We  shall 
instead  examine  in  the  next  section  the  evolution  based  on  an  approximation  to 
the  dispersion  relation  (6). 


III.  Propagation  of  Wavepackets 


With  the  demonstration  that  the  hose  instability  is  convective,  it  is 

meaningful  to  consider  the  spatial  amplification  rate  Im(ajX)  as  a  function  of 

the  frequency  ft.  Shown  in  Fig.  1  is  the  growth  rate  Im(u)X)  as  a  function  of 

real  ft  obtained  from  the  dispersion  relationship  (6)  for  the  spread  mass 

model.  It  is  seen  that  the  presence  of  a  return  current  (f  <  0)  is 
4  5 

destabilizing.  »  The  return  current  yields  a  nonzero  amplification  rate 
at  ft  «  0  when  -f  is  sufficiently  large.  In  fact,  one  can  easily  show  from  (6) 
that  amplifying  waves  are  possible  for  ft-  0  if  -  f>  x/(l  +  X).  Physically, 
the  return  current  has  a  tendency  to  repel  the  beam  axis  from  the  magnetic 
field  axis  and  thus  weakens  the  restoring  force.^'^  Note  from  Fig.  1  that  the 
unstable  amplification  band  is  restricted  to  |ft|  <  0.8.  The  imaginary  part 
of  ( ujX)  is  an  even  function  of  ft  and  the  real  part  of  (a>X)  is  an  odd  function 
of  ft,  ft  being  real. 

The  dispersion  relationship  (7)  is  a  complicated  function  of  ft  which 
defies  analytic  study  of  the  space-time  evolution  of  a  wave  packet.  To  render 
the  problem  analytically  tractable  we  shall  approximate  Im(u>X)  as  a  cubic 
polynomial  in  ft  and  Re(ojX)  as  a  quadratic  polynomial: 

Im(t»  X)  -  a  +  a,ft2  -  a, ft3,  (ft  >  0)  (15) 

a  0  4  j 

ReUaX)  -  -  b2ft2,  (16) 

where  aQ,  a2»  a^,  b2  are  nonnegative  constants.  Curves  of  Im(a>aX)  for  various 
combinations  of  these  constants  are  shown  in  Fig.  2.  Comparing  Fig.  1  with 
Fig.  2,  it  is  seen  that  the  amplification  rate  Im(u)X)  from  Eq.  (7)  can  be 
approximated  quite  well  by  Im(aiX)  of  Eq.  (15).  However,  Eq.  (16)  is  a  poorer 


I 


approximation  to  Re((i>A).  Equations  ( 7 >  — ( 1 0)  show  that  Re(uj\)  reduces  to  a 


sixth-degree  polynomial  in  ft,  for  large  ft.  If  the  beam  profile  were  Gaussian 


rather  than  Bennett,  the  spread  mass  model  would  lead  to  a  fourth  degree 


polynomial. 


To  correct  this  deficiency,  we  shall  bracket  the  real  shape  of  Re(u>A)  by 


allowing  b£  to  range  from  zero  to  infinity.  In  this  way,  we  can  determine 


whether  wave  group  convection  is  sensitive  to  the  model  chosen.  We  also  note 


from  Figs.  1  and  2  that  the  constant  a  in  (15)  may  be  used  to  represent  the 


amount  of  the  return  current,  which  leads  to  a  nonvanishing  value  of 


Im(aiA)  at  ft  -  0. 


If  the  spectrum  S(ft)  of  the  source  is  an  even  function  of  ft,  the  integral 


in  (14)  may  be  carried  over  the  range  0  <  ft  <  »  since  Im(^jA)  is  an  even 


function  and  Re(u>\)  is  an  odd  function  of  ft.  Thus,  instead  of  carrying  out 


the  complicated  integral  (14),  we  are  led  to  consider  the  simplified  integral 


Ra(Z,\J>)  *  /  dftS(ft)e 


-iftZ  -  iAui  (ft)* 


where  oj  (ft)  is  the  approximated  propagation  constant  [cf.  (15)  and  (16)].  If 
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we  further  restrict  consideration  to  the  case  S(ft)  -  1,  substitution  of  (15) 


and  (16)  into  (17)  yields 


a  * 

R  (Z,*)  -  e  °  I(Z,*) 

8 


where 


eo  .  _ 

I(Z,*)  -  /  dft  exp  {-iftZ  +  i},  [(a2+  i  b2)  ft  -  a.jft  ]]  (19) 

o 


is  somewhat  similar  to  an  integral  representation  of  the  Airy  function.  Note 
that  the  effect  of  the  return  current  (aQ)  is  factored  out  in  (18)  and  the 
remaining  integral  I  defined  in  (19)  is  independent  of  aQ.  It  is  readily 
demonstrated  from  (19)  that  I(Z,ijj)  -*•  0  as  Z  +  °»  at  any  fixed  positive  value  of 
Thus,  the  approximate  response  function  Ra(Z,4))  is  consistent  with  our 
earlier  conclusion  that  absolute  instability  is  absent,  regardless  of  the 
return  current. 

The  evolution  of  I (Z , ip)  is  best  visualized  by  setting 

*>  -  VZ  (20) 

with  V  fixed,  but  letting  Z  -►  <*>.  The  result  function  I(Z,VZ)  may  then  be 
interpreted  as  the  evolution  of  disturbances  recorded  by  an  instrument  which 
travels  at  a  normalized  velocity  V  moving  from  the  beam  head  toward  the  beam 
tail.  The  integral  I  is  then  a  function  of  V,  a^,  b2  and  33,  and  the 
asymptotic  response 


lRa(*,*)l  -  e  p*  (21) 

may  be  regarded  as  the  spatial  dependence  recorded  by  the  moving  instrument. 
Again  p  *  p(V)  when  the  constants  aQ,  a2>  aj,  b2  are  fixed. 

We  have  used  the  method  of  steepest  descent  to  evaluate  the  asymptotic 
behavior  (Z  ♦  «>)  of  I(Z,VZ)  as  defined  in  Eq.  (19)  for  various  combinations  of 
aQ,  a2»  aj,  b2  and  V.  The  procedure  is  rather  involved  and  is  outlined  in  the 
Appendix.  Our  findings  are  summarized  below. 

Since  there  are  five  free  parameters,  it  is  convenient  to  Introduce  a 


complex  constant 


=  (a2+  i  b2)/3a3  =  I 1 e 


where  a  =  tan”1  (b2/a2>  ranges  from  zero  to  ir/2  as  b2  varies  from  zero  to 
infinity. 

The  value  of  a  can  be  specified  by  requiring  that  the  value  of  Re(u)X)  or 
its  derivative  at  the  frequency  at  which  maximum  amplification  is 
expected,  match  the  correct  value  from  Eqs.  (7)-(10).  The  value 
of  a  typically  lies  in  the  range  between  15°  and  60°.  Over  this  range,  we 
shall  see  that  the  most  important  features  of  the  results  are  not  sensitive 
to  a;  this  justifies  the  use  of  our  representation  (16). 

Instead  of  using  V,  we  find  it  convenient  to  represent  the  data  in  terms 
of  U,  where 


U  =  SajlG^  V. 


We  normalize  the  spatial  amplification  rate  p  with  respect  to  p^,  the  maximum 
value  of  amplification  rate  expected  from  the  approximation  dispersion 
relation  (15), 


4  a2 


PM  =  Max[Im(uaX)]  -  aQ  +  jj 


Note  also  that  p^  occurs  at  n  ■  where  2a2/3a3,  at  which  the  "group 
velocity"  VM  is  given  by 


3(xHa)|_1 


By  examining  the  asymptotic  response  I(Z,VZ)  as  Z  +  »,  in  a  manner 


teat 
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described  in  Che  Appendix,  the  normalized  spatial  amplification  rate  p/p^  as  a 
function  of  U  for  various  values  of  0^2)  is  shown  in  Fig.  3,  where  we  set  f  ■ 
aQ  *  0.  It  Is  noted  from  this  figure  that  the  maximum  value  of  p  is  indeed 
p^.  This  result  may  be  proved  in  a  rigorous  manner.  Note  also  from  this 
figure  that  there  is  a  minimum  velocity  Umin  for  amplification  to  occur,  which 
is  typically  about  1/3— 1/2  of  the  velocity  UM  for  peak  amplification.  In 
other  words,  an  observer  moving  backward  in  the  beam  at  a  normalized  velocity 
U  which  is  less  than  Um^n  always  sees  the  disturbance  decaying  away  at  late 
times.  It  would  not  be  obvious  from  a  casual  examination  of  the  dispersion 
relation  that  Umin  is  considerably  less  than  Um.  Note  also  that  p  remains 
nonnegative  as  U  -►  «.  This  is  due  to  our  use  of  Lee's  reduced  form  of 
Maxwell's  equations  (appropriate  to  highly  relativistic  paraxial  beams)  to 
derive  Eq.  (7).  In  Lee's  form,  the  field  equations  are  diffusive  when  the 
plasma  is  resistive;  signal  propagation  can  occur  at  infinite  speed  for 
diffusion  equations.  If  the  exact  Maxwell  equations  were  used,  there  would  of 
course  be  no  signal  propagation  at  V  >  c.  This  deficiency  in  Lee's 
formulation  is  of  little  practical  significance.  Shown  in  Figs.  4  and  5  are 
the  values  of  p/pM  for  the  two  cases  f  *  -  0.4  and  -  0.6  [cf.  Fig.  2.]  Again, 
the  minimum  speed  Um£n  of  propagation  for  the  amplifying  waves  is  about  1/2  to 
1/3  of  that  associated  with  maximum  amplification,  for  the  range  of  parameters 
surveyed.  This  mimimum  velocity  is  compared  with  Ujj  *  Sa^ln^l  in  Fig.  6, 
where  %  is  defined  in  Eq.  (25).  Of  course,  Um^n/U^  is  model-dependent,  but 
Umin  <  j  0^  over  the  range  of  parameters  surveyed.  In  fact,  for  low  values 
of  a  (b2)t  is  substantially  less  than  UM  (Fig.  6). 

In  all  cases,  the  value  of  ^min  is  not  particularly  sensitive 
to  the  value  assumed  for  b£>  i.e.  for  a.  Thus,  it  would  appear  that 
our  approximation  (16)  is  adequate,  even  though  there  is  no  well-defined 
way  to  specify  b2> 
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IV.  Remarks 


In  this  paper,  we  have  performed  a  quantitative  study  of  the  evolution  of 

the  convective  hose  instability  based  on  simplified  models  applicable  to  the 

main  body  of  the  beam.  We  find  that  the  instability  is  always  convective  in 

the  beam  frame,  regardless  of  the  return  current.  At  late  time,  the  peak  of 

the  disturbance  convects  backward  in  the  beam  at  a  dimensionless  velocity 

UM  “  I  ao"  1  >  while  amplifying  at  a  rate  P  ■  Im  [  X<d(  SV.)  ] , 

1  'a  -  Gjj 

where  is  the  most  unstable  value  of  Q,  i.e.  the  value  for  which  Imw(jj  )  is 
largest.  However  the  disturbance  also  spreads  out  as  it  convects  backward,  so 
that  the  trailing  edge  of  the  growing  disturbance  moves  at  a  considerably 
slower  velocity  Um£n»  which  we  have  shown  to  be  somewhere  between  1/3  %  and 
1/2  UM.  Along  characteristics  moving  backward  in  the  beam  at  velocities  less 
than  ^min»  the  disturbance  damps  away.  More  precise  specification  of  Uffl^n  is 
not  possible,  within  the  limitations  of  the  simple  model  used.  The  mutual 
influence  of  beam  head  erosion11’14  and  the  convection  of  hose 
wavepackets  remained  to  be  assessed.  As  a  function  of  5  [rather  than  of 
=  X  *  in(s/ro)],  the  instability  grows  fastest  and  convects  slowest  just 
behind  the  pinch  point,  near  the  front  of  the  beam.  It  is  possible  that 
erosion  can  significantly  enhance  instability  growth  when  the  erosion  speed  is 
comparable  to  the  convection  speed  in  this  region,  thus  causing  the  region  of 
the  beam  that  is  most  hose-unstable  to  move  backward  in  5  along  with  the  hose 


wave  packet 
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Appendix 


Asymptotic  Evolution  of  I(Z,4>) 

In  this  appendix,  ve  outline  the  procedure  used  to  evaluate  the  integral 
X(Z s given  by  Eq.  (19)  of  the  main  text,  with  the  understanding  that 
"  VZ,  where  V  is  a  pre-assigned  positive  constant  and  Z  *  ».  The  saddle 
point  method  will  be  used.  Note  that  there  are  two  saddle  points  for  the 
exponents  in  the  integral  of  (19).  Thus,  the  asymptotic  expression  to  (19)  in 
general  contains  three  terms:  the  end  point  contribution  arising  in  the 
neighborhood  of  8  «  0,  and  the  saddle  point  contributions  near  the  two  saddle 
points.  An  inspection  of  (19)  shows  that  the  end  point  contribution  is  less 
significant  than  the  saddle  point  contribution,  and  its  role  will  thus  be 
ignored.  The  location  of  the  saddle  points  determines  which  saddle  point 
contribution  will  dominate  the  Integral  (19),  as  we  shall  show  here. 

With  the  transformation 


n' 


=  <v>1/3  (a  -  nl). 


(A-l) 


where  is  defined  by  Eq.  (22)  of  the  main  text,  it  can  be  shown  that  I( Z , tp ) 
may  be  written  as 


-1/3 

I(Z,f)  -  Const  x  <j,  e  J(Z,^)  (A-2) 


where 


Ej  (Z  ,i|0 


SI 


t  -  iZ)  -  a3<|i 


3 


and 


16 


(A-3) 


i 


% 


=  /  da' 


-  a'3  +  3  a  -2  a' 
s 


In  this  equation. 


n/  =  -  U3*)1/3  a, 


.  .2  _  n  ,2  iZ 

8  '  ^  3(a3t)1/3 


(A-4) 


(A-5) 


and  the  upper  limit  of  the  path  integral  in  (A-3)  approaches  the 

positive  He(a')  axis  (Fig.  7).  From  (A-3),  we  then  see  that  the  saddle  point 

contribution  to  I  depends  just  on  the  complex  quantities  ft.*  and  a  '•  These 

I  8 

latter  quantities,  of  course,  depend  on  a£,  b2,  83,  Z,  <j>,  by  Eqs.  (A-4)  and 
(A-5). 

The  saddle  points  for  J(Z,  \|))  are  located  at  a'  »  agl"  and  a'  "  ag2'  in 
the  complex  a'  plane,  where 


2  2  2 
°.l  “  °s2  “  V  ’ 


(A-6) 


which  is  defined  in  (A-5).  If  either  ag^'  or  ag 2'  *n  a°y  one  t*ie 

shaded  regions  of  Fig.  7,  the  contribution  due  to  that  saddle  point  will  be 
dominant  (over  that  due  to  the  other  saddle  point)  provided  the  contour  for 
(A-3)  may  be  detoured  along  steepest  descent  curves  which  pass  through  that 
saddle  point. 

As  an  example,  consider  the  values  of  a2,  b2»  83,  Z,  such  that  the 
saddle  points  agl'  and  ag2'  are  depicted  as  in  Fig.  8  (left).  The 
contribution  to  J  due  to  a  will  be  subdominant  whereas  that  due 

Si 

to  ag2'  will  be  dominant.  The  steepest  descent  curves  which  pass  through 


and  ns2'  are  also  sketched  In  Fig.  8  (left).  The  initial  point  of 
integration  always  lies  in  the  third  quadrant  of  the  9'  plane,  and  the  path  of 
integration  originating  from  is  also  depicted  in  Fig.  8  (left).  Now, 

consider  domains  A,  B,  C,  D,  all  of  which  lie  in  the  third  quadrant  and  are 
bounded  by  the  steepest  descent  curves  passing  through  ‘  8  right]. 

If  the  initial  point  of  integration  lies  in  A  or  B,  only  the  contribution 
around  will  be  picked  out  by  detouring  the  path  of  integration  as  shown 

in  Fig.  8b.  On  the  other  hand,  if  0^'  lies  within  domains  C,  D,  the  path  of 
integration  is  detoured  so  that  both  saddle  points  and  0g2'  are  passed 

through.  Of  course,  in  the  latter  case,  the  contribution  due  to  &g2' 
would  dominate.  It  is  to  be  noted  that  the  topology  of  the  steepest  descent 
curves  change  as  and  fts2~  vary.  Based  on  a  systematic  examination  of  the 

relative  position  between  and  we  are  able  to  determine  the  relevant 

saddle  point  contributions  to  J.  In  this  way,  the  data  shown  in  Fig.  3-5  of 
the  main  text  is  generated. 
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Figure  1 

The  normalized  spatial  amplification  rate  Im(Xu>)  as  a  function 
of  normalized  frequency  n  for  several  values  of  the  return  current 
fraction  f,  as  given  by  the  spread  mass  model  including  a  self-consistent 
conductivity  model,  Eq.  (7). 
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Figure  2 
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The  approximate  dispersion  relation  Im(Xw  )  =  a  +  a9ft  -  a~ft  ,  Eq.  (15) 

a  o  z  j 


To  match  the  result  of  Eq.  (7)  for  each  of  the  various  values  of  f,  we 


have  chosen  the  following  values  of  aQ,  a2,  a3:  f  -  -  0.6:  (aQ,  a2,  a3) 


(1.1,  21.84,  29.12);  f  -  -0.4:  (0.26,  16.08,  21.44);  f  -  -  0.2:  (0,  12, 


16);  f  -  0:  (0,  8.4,  11.2);  f  -  0.2:  (0,  6.12,  8.16). 


Figure  3 

Value  of  the  amplification  rate  p  as  recorded  by  an  instrument  moving 
backward  in  the  beam  at  scaled  velocity  U,  for  a  case  with  f  »  0.  The 
rate  p  is  normalized  to  the  peak  amplification  rate  fty.  Each  curve 
represents  a  different  choice  of  the  parameter  a  *  arctan 
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Figure  6 


Values  of  Che  scaled  velocity  UM  that  yields  the  maximum  amplification 
rate,  and  of  the  mlmimum  scaled  velocity  ^min  that  yields  any 


amplification,  each  as  a  function  of  the  numerical  parameter  <*.  UM  (a)  Is 


independent  of  f,  while  Umin  (a)  depends  weakly  on  f.  The  ratio 


depends  weakly  on  both  f  and  a« 
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Figure  8 


(a)  Left.  The  contour  r  and  the  steepest  descent  curves  which 


pass  through  saddle  points  fi'  and  n'_. 
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(b)  Right.  Domains  A,  B,  C,  D  in  the  third  quadrant  of  the  plane  in 


which  fl'  lies. 
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